Entanglement detection by Bragg scattering 
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We show how to measure the structural witnesses proposed in [P. Krammer et at, Phys. Rev. 
Lett. 103, 100502 (2009)] for detecting entanglement in a spin chain using photon scattering. The 
procedure, moreover, allows one to measure the two-point correlation function of the spin array. 
This proposal could be performed in existing experimental platforms realizing ion chains in Paul 
traps or atomic arrays in optical lattices. 



Multipartite entanglement plays a crucial role in var- 
ious tasks of quantum information, such as in multi- 
party quantum secret sharing pQ and in the computa- 
tional speed-up in quantum algorithms [2J [3] . This mo- 
tivates the great interest in detecting entanglement in 
many-body systems. An experimentally feasible method 
to achieve this task for spin systems was recently pro- 
posed in Ref. [3] and later tested in a quantum optical 
experiment [5]. Such a method is based on the construc- 
tion of entanglement witnesses related to structure fac- 
tors operators and it is suitable to detect entanglement 
for various families of multipartite entangled states, such 
as the Dicke states, without the need of any a priori 
knowledge on the Hamiltonian describing the physical 
system under consideration. In this paper we discuss 
specifically how this method can be implemented in ion 
chains or atomic arrays in optical lattices. In addition, 
we propose how to measure all elements of the structure 
factor matrix, thereby allowing one to reconstruct the 
two-body density matrix of the spin system. 

We first briefly review the structural witnesses method 
and define the structure factor matrix. Consider a chain 
of N spin 1/2 particles which are ordered in a one- 
dimensional array at the positions r± (i = 1, . . . N). The 
structural witness method is based on the measurement 
of the 3x3 matrix 5(q), which is a function of the trans- 
fered wave vector q and whose elements are given by the 
static structure factors 
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In the above expression S" is the a component of the 
spin operator of particle i at position r\ (a — x, y, z) 
and the average is taken over the initial state of the N 
spins, whose entanglement is to be detected. As shown 
in Ref. [4 , for different values of q and different linear 
combinations of ^"^(q) it is possible to detect different 
families of multipartite entangled states. For example, 
with the choice q = one can construct the witness 
operator Wry 

w D = i- 2 _ (s™(o) + §yy(o) - £«(o)) , (2) 

where 1 denotes the identity operator for the iV-qubit 
system and S aP (q) = J2i<j e^^S? Sf are the struc- 
ture factor operators whose expectation values are given 



in Eq. ([lj). Whenever operator Wd gives a negative 
expectation value on a density matrix of the spin array 
Po, T^{WdPo} < 0, we are guaranteed that the ./V-spin 
state under examination is entangled. The witness in Eq. 
([2]) is suitable for detecting symmetric Dicke states with 
any number of excitations, namely its average value for 
symmetric Dicke states is always negative. 

More generally, the structural witness operators take 
the form 
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where c a are real coefficients with \c a \ < 1, C a (q a ) = 
at(jv-i) (^°(q a ) + S aa {-cf)) and q Q denote three pos- 
sibly different values for the transfered wave vector [5]. 
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FIG. 1: (Color online) An entanglement witness for a spin 
array can be constructed by means of a pump-probe experi- 
ment. The intensity of the scattered light allows one to de- 
termine the expectation value of the witness operator. The 
collection efficiency can be increased when the probe field is 
the mode of a high-finesse resonator [6]. 



The underlying detection method relies on the out- 
come of scattering experiments and is therefore partic- 
ularly suited to detect entanglement in many-body sys- 
tems where local addressing of the individual constituent 
particles is not available. The purpose of this paper is to 
provide a specific implementation scheme, suggesting a 
possible pump-probe experiment one can perform on an 
atomic array or on an ion chain. The setup is sketched in 
Fig. [l] We show that the measurement of the intensity of 
the probe light allows one to extract all quantities needed 
for determining the entanglement witnesses constructed 
through linear superpositions of the quantities in Eq. ([I]) , 



as in the form In addition, it also allows one to deter- 
mine all elements of the structure form matrix, therefore 
measuring the two-point correlation function of the spin 
system and hence giving access to the two-body density 
matrix. 

Let us first denote by \0)j and the two states of 
the spin j = 1, . . . , N at position rj = jdx, where d is the 
interparticle distance and the array is along the i-axis. 
The spins are prepared in the state po and the elements 
ofEq. are evaluated over p . These elements can be 
found by measuring the intensity of a weak probe after 
implementing a dynamics governed by Hamiltonian 
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and = — — Sj). The parameters ao,«i are real- 
valued frequencies, cj> is a phase, q is the difference be- 
tween probe and pump wave vectors, and a the anni- 
hilation operator of the probe field. The dimensionless 
function g(t) gives the temporal form of the excitation, 
which warrants that the dynamics is in the linear re- 
sponse regime, so that the scattered light depends on the 
initial state. Here, g(t) = for t < and max|g(f)| = 1. 
Moreover, the pulse temporal duration is finite and the 
variation is sufficiently smooth to warrant the effective 
dynamics described in Eq. Q. In the following we will 
characterize its characteristic temporal duration by At. 

Since the number of scattered photons at each exper- 
imental run must be small, a high collection efficiency 
of the scattered light is required. This can be achieved 
when the probe is the mode of a high-finesse resonator 
[5] and could be realized in existing experimental plat- 
forms, where atomic arrays [71 [5] or ion crystals [H] have 
been confined and strongly coupled with a cavity field. 
In the following we shall consider that the probe field 
is one mode of a high-finesse ring cavity, which strongly 
couples with an optical dipole transition. 

In order to give a detailed derivation of Hamiltonian 
Q , we first specify the atomic states composing the spin 
transition of the atoms composing the array. We assume 
that the spin of a single atom is composed by radia- 
tively stable states of the hyperfine multiplet of state 
Si/2- This example is taken for convenience, since the 
treatment in the following can be extended to other pairs 
of radiatively stable states. We denote the spin states 
by |0> = |^i /2 ,-l/2) and |1) = |S 1/2 , +1/2), while the 
excited states (which will be needed in order to obtain 
the effective Hamiltonian) are |eo) = \P\/2,~ 1/2) and 
|ei) = |Pi/2,+l/2). For simplicity, we assume that the 
states |0), |1) are degenerate and their energy is set to 
zero, while states |eo), |ei) are both at frequency lj q , so 
that the Hamiltonian for the relevant electronic degrees 
of freedom reads H c \ = Hu3 J2j (l e o)j(eo| + |ei)j(ei|). 
The electronic levels are coupled by cavity field and lasers 
as shown in Fig. [2j The pump lasers, in particular, have 
wave vector kL, frequency ujl, and couple transitions 




FIG. 2: (Color online) A possible excitation scheme, where 
two hyperfine states of the ground state multiplet are the 
spin states. In the scheme here discussed, the realization of 
the witness lies on the capability to adjust the amplitude and 
the relative phase of the lasers with Rabi frequency f2o and 
Qi. 



1 1) j — > |eo)j and |0)j —> \ex)j with maximum Rabi fre- 
quency 51 and fii, respectively. The cavity mode, which 
acts as a probe, is linearly polarized and at frequency oj c 
and wave vector k. We denote by ai and a\ the anni- 
hilation and creation operators of a cavity-mode photon, 
and by and a\ the corresponding operators for the sec- 
ond, degenerate cavity mode at wave vector — k. Both 
modes drive transitions |0)j — > \eo)j and |l)j — > \ei)j 
with vacuum-Rabi frequency g. The Hamiltonian gov- 
erning the coherent dynamics of spins and cavity modes 
reads 

H L = -M c fa|ai + a|a 2 J - ?iA^(|eo)j(e | + |ei)j(ei| 
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and is reported in the reference frame rotating at the 
carrier frequency uil of the laser pulses, with A = lol —ujq 
and S c = uil -w c , respectively. The phase <f> determines 
the relative phase between the fields, since we assumed 
that f2 ,f2i,(? are real- valued. The incoherent dynamics 
due to the radiative decay of the excited state at rate 
7 can be suppressed by driving far-off resonant Raman 
transitions, which corresponds to choosing values of the 
detuning such that |A| 3> g,£lj, \6 c \,K,~f, where n is the 
linewidth of the cavity modes. In this limit the atomic 
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state follows adiabatically the cavity field, namely 

|0)i<eo| = ||0),(0|( ai c ik ^+ a2 e- ik -^) (6) 
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These equations hold at the lowest non- vanishing order in 
the perturbative expansion in 1/|A| and when the width 
of the pulse is such that At 3> 1/|A|. These values are 
substituted in the Heisenberg-Langevin equation for the 
cavity mode, a,£ = — nai + \f2Kag in + [ag, Hi]/{ih), where 
C£,m(i) is the input noise field, satisfying the relations 
Ki n (i)> = 0, (a^ in (t)4,in(*')> = S kt 6(t- t') {j,k= 1,2) 
[10] . Setting a = oi, in the linear response regime the 
solution reads 
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where qi = kL — k, a?o = <?^o/A, cti = g57i/A, and 
f(t) = i J t dre"(' t ^ i ^)( t - r ) £) (r). Note that the detuning 
S' c = 5 C + g 2 /A includes the dynamical Stark shift of the 
spin states [TT] . Equation ^ has been obtained assum- 
ing k \ae |, so that the photons scattered by the atoms 
into the cavity are not reabsorbed but emitted by the 
cavity. The last line of the right-hand side of this equa- 
tion, in particular, corresponds to the dynamics governed 
by the target Hamiltonian, Eq. Q. For the considered 
setup, this is the only relevant term which determines the 
intensity of the probe field, as we will show. 

The intensity of the probe field at the cavity output 
reads J out = (a out (i) t a out (t)), where a out (t) = y/2~Ka(t) + 
ai n (t). It is evaluated using Eq. ([8]), the properties of 
the input operators, and assuming that the cavity modes 
are initially in the vacuum field and can be cast in the 
form J out = 2n\f{t)\ 2 (I + I- m t), where I is the sum of 
the intensities scattered by each atom, while I mt is the in- 
tensity due to the interfence between the scattered waves 
at different sites. In detail, 
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and 
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The coefficients depend on the lasers amplitudes and rel- 
ative phase as follows 
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The possibility to vary these coefficients gives access 
to the general structural witnesses ([3]), as can be easily 
shown using the identity 
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For instance, the element C x (q) is found by setting a y = 
and a x = a, which corresponds to choosing = f2i 
and <j> = 0. 

This setup allows also to measure the off-diagonal com- 
ponents in the matrix Q. For example, by perform- 
ing the measurement for different values of flo, </> and 
on the intensity of the second cavity mode allows one 
to determine the components S xy (q) and S yx (q). The 
other elements involving also the z component, namely 
S z *(q) and S yz (q) (S zx (q) and S xz (q)), are found by 
performing the same measurement procedure after the 
Hadamard transformation^ = ®j(S x +S*)/V2 (Ufj = 
®j{Sj + Sy j) has been made on the qubits of the ini- 
tial state: the density matrix over which one performs 
the trace in Eqs. Q and (10 1 is now p' = U^poU^, 
which corresponds to a Raman pulse composed by clas- 
sical fields [12] . The net effect is hence to transform the 
operators S x into SJ (Sj into SJ) and vice versa, so that, 
for example, Tr{S x Sfp'} = Tr{S|Sf p Q }. In this way also 
the term C z (q) of the general structural witness can be 
accessed, as well as all other elements needed in order to 
reconstruct the two-body density matrix. 

These results show that in the case of Dicke states the 
entanglement witness ^ can be determined by control- 
ling the amplitudes and phases of the driving laser. One 
possibility for extracting the value at q = is when the 
array is along the cavity axis and orthogonal to the laser 
direction, setting dx ■ k = 2im 1 where n is an integer and 
we assumed that the chain is along the x axis. We note 
that localization about the minima is required, so that 
there is no hopping from site to site. The effect of the 
vibrations about the minima will tend in general to lower 
the interferometric contrast, however in the Lamb-Dicke 
regime this detrimental effect is small [T5] . 

In summary, we have presented a specific setup to im- 
plement a recently proposed multipartite entanglement 
detection method [4] in a chain of atoms or ions in an 
optical lattice. The example we provide here shows ex- 
plicitly how the measurement could be performed with 
current-day technology [THS]. These calculations can be 
extended to other setups based on the atomic micro- 
scopes in Ref. 14], or on quantum non-demolition mea- 
surement like the ones discussed in Ref. |15j . We mention 
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that recently an entanglement detection method based on 
the proposal [4] and on a similar protocol |16j has been 
experimentally realized for the motional degrees of free- 
dom of cold atoms in optical lattices [T7] . 

The present entanglement detection scheme is based on 
the measurement of two-point correlation functions, and 
provides, amongst others, a determination of the mag- 
netic susceptibility, which at thermodynamic equilibrium 
and for specific Hamiltonians constitute a macroscopic 
entanglement witness |18j . We conclude by pointing out 
that measurements of higher order correlation functions 



of the scattered field in the proposed system allows one to 
determine higher order correlation functions of the spin 
operators [T5], which are at the basis of witnesses, for 
instance, of GHZ-like multipartite states [20] . 
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